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ABSTRACT 

Jupiter's electron and proton radiation belts are analyzed, with particular 
reference to the sweeping effect of its five inner satellites, located deep within 
its magnetosphere. The characteristics of trapped electrons and protons with a 
magnetic moment of 50 Mev/gauss, considered typical at Jupiter, are calculated. 
The assumption is then made that a particle would be removed from the radiation 
belt if, in its normal motion, it would happen to impact a satellite. The mean 
absorption time before impact is calculated for particles located at the radial 
distance of each of the satellites. This average lifetime was found to be of the 
order of a few days. A characteristic diffusion time near each satellite was 
calculated, assuming violation of the third invariant due to magnetic fluctuations 
associated with fluctuations in the solar wind. This diffusion time was found to 
be long compared with the absorption lifetimes at Europa (L = 9.5), Io (L = 6.0), 
and Amalthea (L = 2.5). Assuming that this process is dominant, these three 
satellites are expected to act as a barrier to the inward diffusion of solar wind 
particles, and fluxes of trapped protons would be reduced by two or three orders 
of magnitude just inside the orbit of Europa; very few could diffuse past Io or 
Amalthea. Due to the ten-degree tilt of the magnetic dipole with respect to the 
rotational axis, however, particles with equatorial pitch angles greater than 69° 
would have significantly longer absorption lifetimes, and a much larger fraction 
could escape absorption by a satellite. We believe that this is the most likely 
explanation for the observed highly-peaked pitch-angle distributions of trapped 
electrons in the region of the decimeter emission, just inside the orbit of 
Amalthea. A flyby mission for which the lowest possible radiation levels are 
desired should be planned so that the spacecraft remains at relatively high 
magnetic latitudes near perijove. 
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1 The 10.4-cm radio brightness distribution measured by 

Berge (1966) with a two-antenna interferometer. The 
peaks of the distribution appear to be about 0 .9 Rj from 
the edge of the disc. The bars in the lower right indicate 


the resolution of the instrument 2 
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magnetosphere. Diameters are not to scale 13 
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4 Sketch used to calculate absorption and diffusion times 

near a satellite with diameter D s 23 

5 Characteristic absorption and diffusion times near each 

satellite. Diffusion dominates near Callisto, but absorption 
dominates at Amalthea, Io, and Europa 26 


6 Sketch of expected off-equatorial proton fluxes at Jupiter, 
assuming diffusion is caused by magnetic fluctuations. 

CRAND is ignored. Fluxes of inward-diffusing protons are 
expected to be reduced by several orders of magnitude at 
Europa (L=9.5) , and none should diffuse past Io or Amalthea. . 36 

7 The fraction F of particles that can diffuse in past Europa 
without being absorbed, plotted as a function of pitch angle. 

The characteristic diffusion lifetime is assumed to be 24 
days. The increase at 69° < a e < 111° is due to the fact 
that during a substantial fraction of these particles' drift 
period , the satellites appear to pass by at latitudes above 

or below their mirror points 44 



INTRODUCTION 


The planet Jupiter is known to have a strong magnetic field and an intense 
belt of trapped electrons. It probably has a proton radiation belt as well, although 
it is not possible to estimate the proton intensities from Earth-based measure- 
ments. The purpose of this paper is to determine some of the probable charac- 
teristics of Jupiter's electron and proton radiation belts, using observations 
where available to determine known quantities, and applying theoretical concepts 
based on our understanding of the Earth's radiation belts to predict characteristics 
which have not yet been measured directly. We give particular consideration to 
the possible effects of Jupiter's five inner satellites, located deep within Jupiter's 
magnetosphere. These satellites exert a sweeping effect, removing trapped 
particles which impact them. Thus it is likely that actual particle fluxes in some 
regions of Jupiter's magnetosphere are much lower than they would be in the 
absence of any satellites. We begin by reviewing our present knowledge of 
Jupiter's radiation belts. 

Trapped electrons. Nonthermal emission in the decimeter wavelength 
region was first discovered by Sloanaker (1959) at 3000 MHz (10 cm). More 
detailed measurements, especially those by Berge (1966) at 10.4 and 21.2 cm 
and by Bransom' (1968) at 75 and 21 cm, clearly established that this decimetric 
radiation is primarily emitted not from the surface of the planet itself, but 
from a region near the equatorial plane, somewhat removed from the surface. 

The 10.4-cm brightness distribution as measured by Berge (1966) is shown in 
Figure 1. Distinguishing characteristics of this radiation component are its non- 
thermal spectrum, the relatively large extent and distinctive shape of the 
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Figure 1. 
feromet 
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emitting region, the relatively high degree of linear polarization (approximately 
22%) but almost complete absence of circular polarization, and the fact that it is 
to some degree beamed in a preferred direction. Extensive reviews of the radio 
emission from Jupiter at both decimeter and decameter wavelengths have been 
published by Roberts and Komesaroff (1965), Warwick (1964, 1967, 1970), and 
Carr and Gulkis (1969). 

Early analyses of this decimetric radiation by Field (1959), Chang and Davis 
(1962), Thorne (1965), and others established quite clearly that it is synchrotron 
radiation from trapped electrons spiraling around the Jovian magnetic field 
lines. Chang and Davis (1962) were able to explain the radiation as being emitted 
from electrons with energies in the 5-75 Mev range and densities of the order of 
10” 3 to 10“ 2 cm -3 , with fields of 0.1 - 1 gauss in the emitting region. In a more 
recent analysis, Warwick (1970) finds that the observed decimeter radiation could 
be produced by 6-Mev electrons with a peak density at 1.8 R } near the equator 
of 1.6 x 10“ 4 cm” 3 and omnidirectional flux of 4.8 x io 6 cm -2 sec -1 , radiating 
in a field of about 2 gauss. The average electron radiation lifetime would be 
about 0.3 years. Luthey and Beard (1972), by comparing the radio emission at 
10.4 and 21-cm wavelengths, concluded that the emission could be produced by 
equatorial electrons with energies of 5-70 Mev and densities of about 1.5 x 10 -3 
L~ 5 ' 4 cm 3 , where L is the radial distance in units of Jovian radii. At the 
Jupiter Radiation Belt Workshop held at the Jet Propulsion Laboratory in 
July, 1971, a workshop model of relativistic electrons was determined ( Divine , 
1971, 1972; Beck , 1972). This model shows a nominal electron belt with a flux 
of 2 x io 7 cm” 2 sec” 1 out to L = 2 and a gradual falloff outside L = 2. An upper 
limit model was a factor of 3 higher in flux. 
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The equatorial pitch-angle distribution of the electrons in the emission region 
must be sharply peaked near a e = 90°, as first pointed out by Chang and Davis 
(1962). Two observational facts lead to this conclusion: (1) the 22% linear polari- 
zation of the decimeter emission, with the maximum of the E-vector lying in the 
magnetic equatorial plane, rather than perpendicular to it; and (2) the beaming of 
the emission in the magnetic equatorial plane. The total beamed intensity is 
observed to drop approximately 11% as observed from a magnetic latitude of 
±13°, as compared with the equatorial intensities ( Roberts and Komesaroff , 1965). 
Electrons with relatively flat helices are required to fit these observations. 

Thorne (1965) analyzed the Roberts and Komesaroff data and concluded that 
a two-component model of the trapped electrons was required to simultaneously 
fit the observed polarization and beaming at 21 cm. One component is relatively 
isotropic, and the other has extremely flat pitch angles. If one assumes an equa- 
torial pitch-angle distribution proportional to sin q a- e , Thorne concluded that the 
first component must have q £ 2 and the second must have q £ 20. He found that 
a distribution proportional to sin 2 a ^ + 2.0 sin^o a ^ was capable of fitting the 
21-cm data. Similar conclusions were reached in the recent study by Clarke (1970). 

No completely satisfactory explanation has been offered for the presence of 
electrons with such sharply-peaked equatorial pitch angles. Chang and Davis 
(1962) pointed out that inward diffusion of particles by violation of the third in- 
variant, while the first and second invariants are conserved, would leak to a 
preponderance of electrons with flat helices. Nakada et al. (1965) calculated 
this effect, but their results do not indicate nearly so much peaking as the obser- 
vations seem to require. In this study we shall suggest an alternate explanation 
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in terms of electrons with flat pitch angles which are much more capable of 
diffusing in past the inner satellites without being absorbed, due to the 10° tilt 
of the magnetic dipole with respect to the rotational axis. 

Chang and Davis (1962) estimated that the time scale for diffusing solar-wind 
electrons into the region of radio emission by violation of the third invariant was 
of the order of a million years, assuming that the diffusion was caused by varia- 
tions in the solar wind strength and changes in the magnetopause position. How- 
ever, the electron lifetime at L — 2 due to synchrotron loss is only about a year. 
Some additional process, therefore, seems to be needed for the electrons besides 
radial diffusion from magnetopause motion, but the nature of this process is 
presently not understood. It may be related to the storm-time processes that 
rapidly energize and move electrons in the Earth's radiation belt. Electrons 
trapped in the Earth's magnetosphere display frequent non-adiabatic time varia- 
tions ( Hess , 1968). They exhibit large temporal fluctuations, short lifetimes, 
large changes due to magnetic storms, and pitch-angle scattering. In his survey 
paper on magnetospherically-trapped particles, Williams (1972) stressed the 
difficulty in obtaining meaningful time averages of electron fluxes, because of 
the great fluctuations in electron intensities observed throughout the trapping 
region. In one well-documented case ( Williams et al. , 1968), relativistic electrons 
appeared to have been injected at L = 3.0 - 3.5, deep within the trapping region, 
during a major magnetic storm. 

Protons in the Earth's radiation belt do not exhibit such rapid temporal 
fluctuations. In general, the electrons show a much more complicated behavior 
than do the protons, and several additional processes act on the electrons that 
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do not operate on the protons. We have little quantitative understanding of these 
additional processes. It would not be surprising, then, to find that Jupiter's 
electrons also behave in unpredictable ways. 

Trapped protons . Nothing is known experimentally about trapped protons at 
Jupiter, because protons radiate synchrotron emission far less efficiently than 
do electrons; they are therefore undetectable by Earth-based measurements. 

We shall, however, review our understanding of the energetic protons in the 
Earth's radiation belt and use this to study Jupiter. 

The terrestrial outer-zone protons first measured by Davis and Williamson 
(1963) from L = 2 to L = 7, with energies between 100 kev and 5 Mev, have a very 
orderly nature and are now quite well understood ( Hess , 1968). Nakada et al. 
(1965) showed that these protons have diffused radially inwards from the mag- 
netopause, conserving the first two adiabatic invariants ^ and J. These protons 
have an exponential spectrum 

j (>E) = k e' E/E ° (1) 

and the characteristic energy of the particles shows a steady, consistent behavior 

E 0 = c L' 3 

for 90° pitch-angle particles at L <6. For other equatorial pitch angles, E q is 
a different but predictable function of L ( Nakada et al. , 1965). The agreement 
between theory and observation shows that radial diffusion without pitch angle 
scattering is the dominant dynamic process for these protons. Nakada and Mead 
(1965) and Tverskoy (1965) calculated equilibrium fluxes of outer-belt protons, 
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assuming a source near the magnetopause, radial diffusion, and atmospheric loss. 
The predicted fluxes agreed rather well with the measurements. Williams (1972) 
concluded that the behavior of both the steady-state and non-steady state outer- 
belt protons could be satisfactorily described by radial diffusion, with a value 
of the diffusion coefficient D equal to 2.4 x 10 9 L 10 R^/ day. 

The energetic protons in the Earth’s inner radiation belt are now also fairly 
well understood. A major source of these protons has long been considered to be 
cosmic ray albedo neutron decay (CRAND). In a recent study Farley and Walt 
(1971) analyzed source and loss processes of inner-belt protons at L < 1.7. 

Using CRAND, energy-loss collisions, and radial diffusion by violation of the 
third invariant, they calculated proton flux intensities, energy spectra, and radial 
distributions. They found good agreement with observations for equator ially- 
trapped protons with first invariant /x between 200 and 3000 Mev/gauss. The 
radial diffusion coefficient required (1.0 x 10~ 8 L 10 R|/day or 6.4 x 10 6 L 10 /x 1 
R 2 /day) was consistent with those derived by a variety of other methods. 

Thus we seem now to have a fairly good quantitative understanding of trapped 
protons in the Earth’s radiation belt. They behave reasonably predictably using 
straightforward concepts of adiabatic theory. The two major sources seem to be 
protons in the plasma sheet whose ultimate source is the solar wind ( Williams , 
1972) and CRAND, with radial diffusion by violation of the third invariant being 
the main mechanism for adjusting their spatial position. 

In this paper we shall assume that protons on Jupiter also behave in pre- 
dictable ways. The known characteristic of Jupiter’s magnetosphere will be 
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examined and some of the characteristics of trapped particles which we might 
expect to find there will be calculated. In particular, the possible effects of 
Jupiter's satellites on its radiation belts will be studied. Our calculations will 
be applied to electrons as well as to protons, although we recognize that a 
number of unknown processes might be operating on electrons. 
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JUPITER’S MAGNETIC FIELD AND MAGNETOSPHERE 

30 

The strength of Jupiter's magnetic dipole moment is approximately 4 x 10 
gauss-cm 3 , or 10 gauss-R 3 (Warwick, 1970). Although several independent 
arguments lead to a number of this magnitude, perhaps the most convincing one 

O 

relates to the explanation of the extremely sharp upper-frequency cutoff of the 
Io-associated decametric emission at 39.5 MHz. If we assume that this cutoff 
frequency is equal to the electron gyrofrequency at the base of a dipole field 
line passing through Io (L = 5.95), we get a surface field of 14.1 gauss at a 
latitude of about 65°, or somewhat less than 10 gauss at the equator. For the 
calculations in this paper, we shall assume a field strength in the magnetic equa- 
torial plane equal to 10/L 3 gauss. At corresponding L-values, this is about 30 
times greater than the Earth's field. 

The Jovian dipole is tilted by about 10° with respect to the rotational axis. 
Measurements of the decimetric radiation by Roberts and Komesaroff (1965) show 
that the plane of polarization rocks back and forth ±10° relative to the rotational 
equator as the planet rotates. This indicates that the magnetic equator is in- 
clined by about 10° with respect to the rotational equator, and the magnetosphere 
and its trapped radiation exhibit a diurnal wobble similar to the Earth's. Careful 
measurements by Berge (1972) of the position of the decimetric radiation with 
respect to the visible planet indicate that the center of the dipole is offset, 
vertically and horizontally, by less than 0.1 Rj. 

The high magnetic field strength at Jupiter leads to an extremely large 
magnetosphere. From considerations of pressure balance, the distance to the 
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magnetopause is proportional to (B 2 /m n v 2 ) 1/6 , where n is the solar wind 
particle density and v its velocity. If the solar wind extends to Jupiter, the flux 
n v must be inversely proportional to the square of the distance from the sun. 

If we assume that the velocity is constant, the particle density at Jupiter is down 
by a factor of (1/5.2) from that at Earth. Thus the ratio of magnetopause distances 
in units of planetary radii will be 

rjupiter = (5.2 x 

For the calculations in this paper we shall assume a Jovian magnetopause at 50 Rj 
in the solar direction. The magnetospheric cross section as viewed from the 
sun or Earth would be about 200 Rj in diameter. Jupiter's magnetosphere at 
opposition would therefore have an apparent diameter in the sky of about 1.3° 
and an apparent cross section area about seven times as great as the sun or moon. 

The magnetospheric plasma and associated trapped radiation co-rotate with 
the planet out to about 50 Rj ( Carr and Gulkis , 1969; Brice and Ioannidis , 1970). 
This conclusion is reached from an analysis of the formation of the Earth's 
plasmapause, which marks the boundary of the co-rotating portion of the mag- 
netosphere ( Brice , 1967). The plasmapause occurs where the convective electric 
field is about equal tothe co-rotation electric field. Whereas at Earth they are 
equal at about 5 R E , on Jupiter the co- rotation electric field is everywhere orders 
of magnitude greater than the convective field ( Brice and Ioannidis , 1970). We 
therefore expect to find co-rotation of trapped particles all the way out to Jupiter's 
magnetopause. 


JL) 

0.31/ 


1/3 


= 5.5 


(3) 
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CHARACTERISTICS OF JUPITER'S INNER SATELLITES 


Jupiter has twelve known satellites. Of these, the five inner ones might be 
expected to have a significant effect on its trapped radiation. These consist of 
the four Galilean satellites, large enough to be seen even with small telescopes, 
plus the small satellite JV, sometimes referred to as Amalthea, nearest to 
Jupiter. The remaining seven satellites are very small, much further away from 
Jupiter, and are in highly inclined and elliptical orbits. 

The basic characteristics of the five inner satellites are given in Table 1. 

All five are in almost precise circular orbits with zero inclination. L-values 
at the position of each satellite are calculated, using a Jovian equatorial radius 
Rj of 70,850 km ( Dollfus , 1970a). Due to the large oblateness of Jupiter, these 
L-values would be about 2.2% larger if a mean radius (69,320 km) were used 
instead of the equatorial radius. Figure 2 shows the position of these satellites 
with respect to Jupiter and its magnetosphere. 

Also calculated is the co-rotation period of each satellite. This is the 
apparent (retrograde) period in a frame of reference rotating with Jupiter's 
decametric rotation period, 9 h 55 m 29.73 s ( Donivan and Carr , 1969). The decametric 
period is believed to be associated with Jupiter's solid surface and internal 
magnetic field; it is used to define Jupiter's System HI longitude system. The 
satellites' Systemlll co-rotation period will be important in a later section where 
we analyze the effect of satellite motion relative to trapped particles co-rotating 
with the decametric period. 
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TABLE 1 


Characteristics of Jupiter’s five inner satellites. The radial distance, sidereal 
period, inclination and orbital eccentricity are taken from Blanco and McCuskey 
(1961). L is the radial distance in units of one Jupiter equatorial radius, Rj = 
70,850 km ( Dollfus , 1970a). The co-rotation period r cr is the apparent (retro- 
grade) period of the satellite in a frame of reference rotating with the decametric 
rotation period, 0.413539 days ( Donivan and Carr , 1969). The satellite diameters 
are taken from Dollfus (1970b). 

SATELLITE 



JV 

JI 

JII 

jni 

JIV 


Amalthea 

Io 

Europa 

Ganymede 

Callisto 

Radial 

distance, kms 

180,500 

421,600 

670,800 

1,070,000 

1,882,000 

L 

2.55 

5.95 

9.47 

15.10 

26.56 

Inclination 

0°24' 

0°2’ 

0°28’ 

0°11’ 

0°15’ 

Eccentricity 

0.003 

0.0000 

0.0003 

0.0015 

0.0075 

Sidereal 
period, days 

0.498 

1.769 

3.551 

7.155 

16.689 

Co-rotation 
period r cr , days 

2.434 

0.540 

0.468 

0.439 

0.424 

Diameter 
D s , km 

200 

3500 

3100 

5550 

5000 
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L3 


The diameter of Io is about the same as the Earth's moon, Europa is slightly 
smaller, and Ganymede and Callisto are about 50% larger. Amalthea is very 
much smaller, but due to its innermost position it can significantly influence the 
radiation belts. 
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TRAPPED-P ARTICLE MOTION AT JUPITER 

We now examine the behavior of trapped particles in Jupiter's magnetosphere. 
We assume that particles co-rotate with the planet and that the equatorial magnetic 
field is 10/L 3 gauss. Typical characteristics of near-equatorial protons and 
electrons near the surface (L = 1) and near each of the five inner satellites are 
shown in Table 2. We assume conservation of the magnetic moment fj.. Thus 
particle energies increase with decreasing distance. 

Nonrelativistically, the kinetic energy of a particle is given by E = i uB m , where 
B m is the field magnitude at the mirror point. However, the electrons in the 
inner regions will be relativistic, and the relativistic form of the first invariant 
must be used: 

M = p*/2 m 0 B = p 2 /2 m 0 B m (4) 

where p is the relativistic momentum and m Q is the rest mass. The kinetic 
energy is then given by 

E = (p 2 c 2 + m 2 c 4 ) 1/2 - m 0 c 2 . fi) 

In the extreme relativistic limit this reduces to 

E=(2tn 0 c 2 M B m ) 1 / 2 -m 0 c 2 ( 6 ) 

with m 0 c 2 = 0.511 Mev for electrons. Thus, very roughly, E 2 = m B m for 
relativistic electrons, with energies in Mev. 

We have chosen a value of 50 Mev/gauss for the magnetic moment of both 
protons and electrons as characteristic of typical particles which might become 
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TABLE 2 


Characteristics of near-equatorial protons and electrons near the surface of 
Jupiter (L = 1) and at the distance of each of the five inner satellites, assuming 
a magnetic moment M = 50 Mev/gauss and an equatorial surface field at Jupiter 
of 10 gauss. Eis the kinetic energy, R c is the cyclotron radius, T b is the bounce 
period, and r d is the drift period. 



L = 1 

Amalthea 

Io 

Europa 

Ganymede 

Callisto 

L 

1.00 

2.55 

5.95 

9.47 

15.10 

26.56 

B, gauss 

10.00 

0.6048 

0.0475 

0.0118 

0.0029 

0.00053 

Protons 

E, Mev 

410 

29.8 

2.37 

0.589 

0.145 

0.027 

R , km 

C 

3.2 

13.1 

46.9 

94.1 

190 

442 

r b , mins 

0.016 

0.122 

0.991 

3.16 

10.1 

41.6 

r d , days 

0.35 

1.63 

8.65 

21.8 

55.6 

172 

Electrons 

E, Mev 

22.1 

5.07 

1.128 

0.418 

0.129 

0.026 

R , km 

C 7 

0.075 

0.31 

1.09 

2.20 

4.42 

10.3 

T b , mins 

0.012 

0.030 

0.074 

0.13 

0.30 

1.02 

r d , days 

10.8 

17.3 

27.7 

39.7 

69.6 

181 
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trapped in Jupiter's outer magnetosphere. As indicated in Table 2, maximum 
equatorial proton energies for this value of m are about 400 Mev. Electron 
energies in the heart of the synchrotron emission region are about 10 Mev, con- 
sistent with earlier estimates. Brice (1972) estimates that protons inside Jupiter's 
bow shock would have a magnetic moment of about 100 Mev/gauss, and electrons 
about 20 Mev/gauss. 

Note that if the first invariant is conserved, the kinetic energy of relativistic 
equatorial electrons varies as L 1-5 . whereas the energy of nonrelativistic protons 

“3 

varies as the usual L . Protons are therefore much more energetic than elec- 
trons with the same magnetic moment in the inner regions, where for electrons 

aJB » m Q c 2 . 

Shown also in Table 2 are the cyclotron radius R c , the bounce period r b , and 
the drift period r d inthe co-rotating frame for near-equatorial particles, as 
calculated from the following formulas ( Hess , 1968), applicable to particles of 
any energy: 


R c - P c/e B 

(7) 

r b = 3 L Rj/v 

(8) 

e B Q Rj/3 L m 0 v 2 c 7 

(9) 


where v is the particle velocity, and B 0 is the surface equatorial field (10 gauss). 
The drift period scales as B 0 R 2 , where R is the planetary radius, and therefore 
is several orders of magnitude larger on Jupiter than on Earth. In Jupiter's 
field, the drift period for near -equatorial nonrelativistic protons or electrons 
in the co-rotating frame simplifies to r d = 122/LE days, with E in Mev. The 


17 



drift period for highly relativistic particles is a factor of two greater. For non- 
equatorial particles r d is multiplied by an additional factor F(\), ranging from 
1.0 to 1.5, where X. is the mirror latitude. The bounce period for non-equatorial 
particles is also multiplied by a factor ranging from 1.0 to 1.9, depending on the 
mirror latitude. Thus the values given in the table are valid within a factor of 
two for non-equatorial particles as well. 

Nonrelativistically, R c a , r b a M 'i /2 , and r d <x M -1 .In the extreme 
relativistic limit R c a M 1/2 , r h is constant, and T d a M~ 1/2 . Thus the values 
given in Table 2 are not strongly dependent on the specific value we have chosen 
for the magnetic moment. 

A comparison of Tables 1 and 2 leads to the following conclusions: 

1. Cyclotron radii are small compared with satellite diameters. Thus we 
need consider only the motion of the guiding center as we study the interaction 
of trapped particles with Jupiter’s satellites. 

2. Drift periods in the co-rotating frame are very long on Jupiter, and 
co-rotation dominates. Jupiter's energetic trapped particles thus behave very 
differently from those at Earth, where drift periods are typically short compared 
with the 24-hour rotation period. 

We can now form a picture of particle motion in the vicinity of the five inner 
satellites. Since the relative motions are fairly complex, one can easily reach 
erroneous conclusions unless one thinks in terms of the most appropriate co- 
ordinate system. We have found that the best system to visualize is a frame of 
reference as shown in Figure 3. This system is centered on the magnetic dipole, 
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with its z-axis aligned with the dipole, and co-rotatingwith the planet. The axes 
are fixed with respect to the solid planet and define the jovimagnetic dipole 
coordinate system. In this system the low-energy plasma appears to be stationary, 
and high-energy trapped particles drift in longitude with periods of the order of 
weeks. Jupiter is stationary, and its satellites appear to be in retrograde orbits 
with an inclination of 10° with respect to the equator and periods given by the 
co-rotation period in Table 1, i.e., about half a day except for Amalthea. The 
satellite orbits, although tilted, are coplanar, with their common nodes defined 
by the intersection of the jovimagnetic equator with the rotational equator. These 
nodes do not precess. The satellites appear to differ only in their radial distance 
and their apparent rate of retrograde motion. Due to the offset of the dipole, the 
satellite orbits, although circular, appear to be offset by an equal amount; perijove 
is perhaps 0.2 Rj closer than apojove. Perijove, however, is always at the 
same jovimagnetic longitude for all satellites. The plane in which the satellites 
appear to revolve is perhaps as much as 0.1 Rj above or below the jovimagnetic 
equator, in addition to being tilted by 10°. 

In the next section we shall use a simplified version of this picture, where 
we ignore the effects of dipole inclination, dipole offset, quadrupole and higher- 
order terms, and distortion of Jupiter’s magnetosphere topology by the solar 
wind. In later sections we will examine the consequences of removing some of 
these simplifying assumptions. 
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Figure 3. Co-rotating jovimagnetic coordinate system. Jupiter and the thermal plasma are stationary in this system. The inner 
satellites appear to revolve in a retrograde direction with an inclination of 10° and periods of about 10-13 hours 
(60 hours for Amalthea). 



SWEEPING EFFECT OF JUPITER'S SATELLITES 


The possibility that a planetary satellite might sweep out a region of a radia- 
tion belt was first considered by Singer (1962). At the time it was thought that 
Mars had a magnetic field and a radiation belt similar to the Earth's. Singer 
calculated the sweeping effect of the tiny satellite Phobos, located at a martian 
L = 2.8. In this section we calculate the rate at which near -equatorial trapped 
particles are removed by Jupiter's inner satellites. To simplify the analysis, 
we make the following assumptions: 

1. The dipole is centered and aligned with the rotational axis. 

2. The field is purely dipolar — i.e., no quadrupole or higher-order effects 
and no external distortions. 

3. The particles' guiding-center motion is not affected by the nearby presence 
of a satellite. Thus we exclude the possibility that a particle might be deflected 

by the electric or magnetic field surrounding a satellite. If in its normal motion 
a particle would impact a satellite, we assume that the particle will be removed 
from the radiation belt. This is the most crucial assumption, and if it is not 
true, the major conclusions of this study would be substantially altered. 

As in the previous section, we assume that all trapped particles co-rotate 
with the planet. In our co-rotating reference frame the satellites will appear to 
revolve in the equatorial plane with retrograde co-rotation periods t cr given by 
Table 1. As the particles bounce back and forth through the equator, they have 
a finite probability of impacting a satellite. In the analysis we distinguish between 
two cases, depending upon the longitudinal distance d that a satellite appears 
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to move with respect to the trapped particles during one half bounce period. 

This distance is given by 

d = tt L R t t /t (10) 

J b' cr 

If this distance is greater than a satellite diameter D s (case 1), the absorption 
is a random process; in one apparent revolution of the satellite only some frac- 
tion of the particles at the position of the satellite will be absorbed. On the 
other hand, if d is less than a satellite diameter (case 2), the satellite will act 
like a snowplow and remove all particles in its path. 

Case 1; Random absorption (d > D s ). This case was considered by Mead 
(1972). We consider a narrow ring of monoenergetic trapped particles whose 
radial width in the equatorial plane is exactly the diameter of the satellite D 

S 

(Figure 4). For the moment we assume no diffusion of particles into or out of 
the ring. In one half bounce period each particle will pass through the equatorial 
plane once. Those which happen to impact the satellite will be removed. One 
half bounce period later the satellite will appear to have moved in longitude by 
at least one satellite diameter; therefore, it will absorb a whole new group 
of particles. Many of the particles in between will be missed, but these will 
have a finite probability of being absorbed the next time around. 

In one half bounce period, r h /2, the average probability P that any one 
particle in the ring will be absorbed is therefore given by the ratio of the 
absorption cross section of the satellite (tt D 2 s /A) to the area of the ring (2 tt L 
Rj D s ); i.e., 

P = D /8 L R t ( n ) 

s J 
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Figure 4. Sketch used to calculate absorption and diffusion times 
near a satellite with diameter D s . 
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The fraction lost per unit time is then given by 2P/r b = D s /4 r fa L R J} and over 
many co-rotation periods the mean absorption lifetime r abs before impacting 
the satellite is given by the inverse of this quantity, i.e., 

r abs = V 2 P = 4 r b L Rj/D s (Case 1) (12a) 

Since r a for nonrelativistic' particles, the Case 1 absorption lifetime 

varies as the inverse square root of the magnetic moment. 

Case 2: Snowplow absorption (d < ). If the apparent drift of the satellite 

in one half bounce period is less than a satellite diameter, all particles in the 
center of the ring (Figure 4) will be absorbed in one corotation period. Although 
some near the edge of the ring might slip past, the mean lifetime r abs will be 
approximately equal to the co-rotation period: 

T abs * T cr ( CaSe 2 ) (12b) 

The Case 2 absorption lifetime is independent of the magnetic moment of the 
particle. Except for a factor of 7t/4, the Case 1 and 2 absorption lifetimes are 
equal for d = D s . 

In Table 3 values of d are given for electrons and protons with magnetic 
moment p. = 50 Mev/gauss at the position of each of the satellites. It is clear 
that for protons, Case 1 applies at Europa, Ganymede, and Callisto, and Case 2 
applies elsewhere. For electrons, which bounce through the equator much more 
frequently due to their greater velocity, Case 1 applies at Callisto, but Case 2 
applies elsewhere. The values of the mean lifetime before absorption t are 
given in Table 3 and plotted in Figure 5 for electrons and protons, with the 
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TABLE 3 

Comparison of absorption and diffusion rates for particles with ^ = 50 Mev/ 
gauss. The parameter d is the longitudinal distance each satellite appears to 
move with respect to the co-rotating system during one half bounce period; it de- 
termines whether Case 1 or Case 2 applies, r is the mean lifetime before 
absorption and r d . ff is a characteristic diffusion time near each satellite. 
A r is the inward distance that an equatorial particle attached to afield line would 
move during a sudden commencement which moved the boundary in from 50 to 
40 Rj. F is the fraction of the particles which survive absorption in diffusing 
past each satellite. ' 
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Satellite 



Amalthea 

Io 

Europa 

Ganymede 

Callisto 

d, protons 

20 km 

1700 km 

9900 km 

54,000 km 

403,000 km 

d, electrons 

5 km 

125 km 

420 km 

1600 km 

9900 km 

T abs • protons 

2.43 d 

0.54 d 

1.90 d 

5.43 d 

43.5 d 

r . , electrons 

abs 7 

2.43 d 

0.54 d 

0.47 d 

0.44 d 

1.07 d 

Magnetic diffusion 

T di f f 

137 yr 

8.7 yr 

24 d 

0.7 d 

0.002 d 

Ar 

10 km 

308 km 

2090 km 

14,800 km 

165,000 km 

F, protons 

10" 125 

10' 66 

0.0016 

0.788 

0.999 

F, electrons 

10" 12S 

10-66 

1.2E-6 

0.159 

0.996 

Electric diffusion 

T 

208 d 

393 d 

19 d 

3.7 d 

0.1 d 

di f f 

F, protons 

10 8 

10" 23 

0.0036 

0.370 

0.995 

F, electrons 

10 " 8 

10" 23 

6.0E-6 

0.006 

0.838 
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CHARACTERISTIC Tl 


AMALTHEA 10 EUROPA GANYMEDE 


CALLISTO 


0.01 d 



L (Rj) 


Figure 5. Characteristic absorption and diffusion times near each satellite. Diffusion dominates 
near Callisto, but absorption dominates ot Amalthea, Io, and Europe. 




appropriate Case 1 or 2 formula being used. These lifetimes are of the order of 
a few days, except for protons near Callisto. The values are rather insensitive 
to the magnetic moment of the particles. 

There will be a group of particles mirroring at latitudes less than 10° (equa- 
torial pitch angles greater than 69°) that will not strike the satellites as often 
because of the 10° tilt between the magnetic axis and the rotational axis (see 
Figure 3). The mean absorption lifetime of these particles will be considerably 
longer than for the particles which mirror at higher latitudes. We will consider 
these particles in a later section. 

Whether or not particles can remain in the radiation belts, despite their being 
continually swept out by the satellites, depends on how rapidly radial diffusion 
moves them past the orbit of the satellites, and, in the case of electrons, whether 
additional source mechanisms are present. In the next section we attempt to 
compare rates of radial diffusion with rates of absorption by the satellites. 
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RADIAL DIFFUSION RATES AT JUPITER 


In an earlier section we argued that radial diffusion by violation of the third 
invariant was likely to be the main mechanism for moving solar-wind protons 
into the inner regions of Jupiter's magnetosphere. Their energy increases with 
decreasing distance from the planet if the first two invariants are conserved. 
Trapped electrons would also undergo radial diffusion. However, a number of 
additional processes, such as pitch-angle diffusion and local acceleration 
mechanisms, are likely to affect the electrons. The diffusion analysis in this 
section is therefore expected to apply primarily to protons, which we believe 
should behave in predictable ways. The behavior of electrons, as at Earth, is 
likely to be much more complex. 

The most difficult part of any trapped-particle diffusion analysis is to 
estimate the value of the diffusion coefficient D, which governs the rate of 
diffusion in different regions of the magnetosphere. At Earth both theoretical 
and experimental techniques have been used to estimate D • In recent reviews, 

Walt (1971a, b) has summarized the present status of radial diffusion studies. 
Almost all the studies indicated a diffusion coefficient increasing extremely 
rapidly with L. The theoretical studies led to a coefficient proportional to L 10 
or L 6 , depending upon whether magnetic disturbances or electric potential disturb- 
ances were considered to be the dominant diffusing force. The limited experi- 
mental studies generally could not determine the exact dependence upon L, except 
for one study where D was found to vary as L 10 ± 1 ( Newkirk and Walt , 1968). 

There was in general a scatter of several orders of magnitude at all L-values 
in the estimated value of D. This was understandable, considering the variety 
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of ways in which the values were obtained, and the fact that both proton and 
electron studies were considered together. In general, the various analyses 
have not considered the possible dependence of D upon the particle energy or 
species. One exception is the recent theoretical analysis by Cornwall (1972) 
of the radial diffusion of ionized helium and protons, where the dependence of 
the electrostatic radial diffusion coefficient upon the charge Z, mass number A, 
and energy was specifically considered. 

At Jupiter the problem of estimating D is even more difficult. We have only 
indirect estimates of electron energies and intensities, and no experimental in- 
formation on protons. We have almost no way of obtaining quantitative estimates 
of fluctuating magnetic and electric fields, except to assume that solar wind 
fluctuations at Jupiter are similar to those at Earth. In this section we take a 
very simplified approach to the determination of diffusion rates at Jupiter. We 
treat only equatorial particles, and we assume that D is independent of particle 
energy. We attempt to estimate a magnetic diffusion coefficient, using an approach 
similar to that taken by Nakada and Mead (1965) in their study of diffusion of 
protons in the Earth's radiation belt. We then compare typical magnetic diffusion 
times in the vicinity of each satellite with the absorption times calculated in the 
previous section, to see whether particles can diffuse past the satellite before 
they are absorbed. And finally, we examine the possible effects of an L 6 -type 
diffusion, as might be produced by electric field fluctuations. 

Magnetic diffusion. Here we estimate the rate of radial diffusion at Jupiter 
due to magnetic fluctuations associated with fluctuations in the solar wind. We 
restrict our analysis to equatorially-trapped particles. As in Nakada and Mead 
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(1965), we assume violation of the third invariant associated with magnetic dis- 
turbances caused by sudden increases in the solar wind intensity, followed by a 
gradual relaxation of the magnetosphere, during which the third invariant is 
conserved. 

We assume that the distorted field strength in the equatorial plane is given by 

„ M 3 1 a 2 , / -« o\ 

B = — + — + — r cos cp 

r 3 r 3 r 4 

b b 

where M is the surface equatorial field strength (10 Gauss), r is the distance in 
jovian radii, r b is the distance to the magnetopause in the solar direction, 0 is 
the solar magnetic longitude (<fi = 0 at magnetic local noon), and from Mead’s 
(1964) model a x = 0.816 M and a 2 = 0.673 M. 

Near Jupiter, the quiet -time contribution to the total field from the magneto- 
pause currents is given by a t /r 3 , where r b is the quiet-time boundary distance, 
which we assume to be 50 . If a sudden increase in the solar wind intensity 

moves the boundary in to a new position r c and the magnetosphere subsequently 
relaxes slowly, the mean square radial displacement of a particle at a distance 
r is given by 


((A r) 2 > 




(14) 


( Nakada and Mead , 1965, Equation 3). 

We now assume that on the average there is one disturbance per (Earth) day 
which doubles the contribution of the external field near Jupiter, i.e., 
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(15) 


V r c = 2 V r b 

This is equivalent to a doubling of the quantity n 1/2 v at Jupiter, where n is the 
solar wind density and v is its velocity. Then the resulting average diffusion 
coefficient D is given by 




25 

196 




(2 4/3 - l) 2 


(16) 


= 0. 13 r 10 /r® day -1 
= 3.4 x 10“ 15 L 10 Rj/day 
: = 2.0 x 10“ 10 L 10 km 2 / s ec 


taking r b = 50 Rj. 

In Nakada and Mead (1965) the diffusion coefficient was determined by count- 
ing the number of sudden commencements and sudden impulses of varying sizes 
observed per year. This led to a value of D equal to 0.0155 r /r® day , or 
1.55 x 10“ 10 L 10 R 2 /day, taking r = 10 R . Thus if an analysis similar to the 
present one had been used in that paper (one sudden doubling of the solar wind 
velocity per day), the resulting diffusion coefficient would have been about ten 
times as large, and roughly equal to the value quoted by Williams (1972) as that 
value which could explain the behavior of outer-belt protons. Nakada and Mead 
(1965) indicated that an increase in their diffusion coefficient by a factor of 8 
would have yielded better agreement with experimentally-measured outer-belt 
proton fluxes, and Falthammar (1966) suggested that consideration of the effect 
of a large number of small, nearly continuous field fluctuations could give the 
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required increase in D. In the present study, therefore, although we have probably 
overestimated the frequency of large disturbances at Jupiter, consideration of a 
spectrum of magnetic disturbances would very likely lead to a value of D similar 
to the value used here. 

We now wish to compare these diffusion rates with the loss rates due to ab- 
sorption by the satellites. The simplest way is to calculate a characteristic 
diffusion time in the vicinity of each of the satellites and compare this with the 
mean lifetime before impact calculated in the previous section. 

A characteristic diffusion time can be defined as follows. The one-dimensional 
Fokker-Planck equation for radial diffusion of trapped particles by violation of 
the third invariant is given by 


3 n _ 1 3 

3T” 2 3~r 



r 2 3 r 





+ S (r, t) 


(17) 


( Birmingham et al. , 1967, Equation 2.4), where n(r, t) is the ensemble-averaged 
number of particles in the range r to r + dr at time t, and S is the source and/or 
sink function describing injection and losses. Let us assume that in the immediate 
vicinity of the satellite S = 0 and that D is independent of r. The equation is then 


3 n 
3T 






3 n 2 n\ 

T7 + ~) 


(18) 
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We now assume that at time t = 0 there are n 0 particles all located at a 
distance r = r 0 . So long as the width of the evolving distribution function is 
narrow compared with the distance r Q ,Bn/Br » 2n/r, and thus 


3 n _ n 3 2 n 


( 19 ) 


which is in the form of the familiar diffusion equation. The solution satisfying 
the initial boundary condition is 


n (r, t) = ° exp (- (r - r ) 2 /4 D t) (20) 

(4 77 D t) 1/2 0 

which can be verified by substitution into Equation 19. Equation 20 describes the 
early evolution of a delta function of n Q particles at r = r Q , t = 0 into a gaussian 
of increasing width. 

We now suppose that the n 0 particles are just at the outer edge of the ring 
of width D s (see Figure 4). We then ask how long it will take for n(r, t) to 
develop into a gaussian whose value at the inner edge of the ring (r = r Q - D ) 
is just 1/e of its maximum value at the outer edge. At the inner edge 

( r - r 0 ) 2 /4 D t = 1 (21) 

and therefore 

7 diff =D'/4D (22) 

The calculated values of this characteristic diffusion time at each of the five 
inner satellites are given in Table 3 and plotted in Figure 5. It is evident that 
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the diffusion times are long compared with the absorption times at Amalthea, Io, 
and Europa. The extreme variation from Amalthea to Callisto is due, of course, 
to the L 10 dependence of the diffusion coefficient. 

The fraction of particles that can diffuse in past a satellite without being 
absorbed can be calculated by obtaining a steady-state solution to the Fokker- 
Planck equation in the region of satellite absorption. We assume that the rate 
of absorption is constant throughout the ring within which the satellite moves 
(see Figure 4) and can be characterized by an average absorption time r abs . 
Then the sink term S(r, t) = n/r abs and the Fokker-Planck equation is 


2 n _ p 2 2 n 



T abs 


(23) 


where, as above, D is assumed to be constant throughout the ring and dn/d r > > 
2n/r. It is convenient to shift the origin by setting r' = r - r Q + D s , whereby r’ = 0 
at the inner edge of the ring. Equation 23 is unchanged and its general solution 
is 


n (r') = A e r /a + B e -r /s 


(24) 


where a =(Dr abs ) 1/2 is a characteristic diffusion length which depends both on 
the diffusion coefficient and the average absorption time. As boundary conditions, 
we let n be arbitrary at r* = D s (the outer edge of the ring) but set n =0atr'=0 
(the inner edge). This is equivalent to placing a sink at the inner edge of the 
ring, which then prevents the diffusion of particles back into the region of 
absorption. For these boundary conditions B = - A and 


n (r') = A (e r ' /a - e" r ' /a ) = 2 A sinh (r'/a) 


(25) 
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The radial flow of particles is given by 


J (r') = - D — = - A_l(e r ' /a + e" r ' /a ) (26) 

3 r ' a 

where a negative value of J represents a net inward flow of particles. 

The fraction F of particles that can diffuse in past a satellite without being 
absorbed is given by the ratio of the inward flow at the inner and outer edges: 

F _ J ( r ' - 0) = 2 (27) 

J( r ' = D S ) e D ./* + e -D,/" 

= 1/cosh (D s /a) = 1/cosh (2 ^ r di{f / T abs ) 

since r diff /r abs = (D s /2a) 2 . Thus, if the characteristic diffusion length a is 
short compared to a satellite diameter, or equivalently, if the absorption time 
T abs * s short compared to the diffusion time, only a small fraction of the 
particles will survive. This fraction is given in Table 3 for particles near each 
of the inner satellites. 

It is clear from the values of T abs , r dif f , and F in Table 3 and Figure 5 
that Callisto will have essentially no effect on trapped particles. However, if our 
value of D is reasonably accurate, given the assumptions of our study, Ganymede 
will exert a relatively small effect, Europa will reduce the flux of inward- 
diffusing particles by several orders of magnitude, and Io and Amalthea will act 
as complete barriers to diffusion. As indicated earlier, we expect that the 
diffusion analysis would be primarily applicable to trapped protons. A sketch 
of the expected proton fluxes at Jupiter, therefore, is shown in Figure 6. An 
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FLUX (Arbitrary Units) 



Figure 6. Sketch of expected off-equatoria! proton fluxes at Jupiter, assuming diffusion is caused by magnetic fluctuations. 
CRAND IS- ignored. Fluxes of inward-diffusing protons are expected to be reduced by several orders of magnitude at 
Europa (L-9.5), and none should diffuse past Io or Amalthea. 



implicit assumption, of course, is that all trapped protons diffuse in from the 
outer magnetosphere and that other sources,, such as cosmic ray albedo neutron 
decay (CRAND), can be ignored. 

We might ask whether particles can be moved radially past a satellite in one 
large jump by a single large sudden commencement. In their study of radial 
diffusion, Davis and Chang (1962) derived an equation for the intersection of a 
distorted line of force with the equatorial plane in a field given by Equation 13, 
where the solar wind moves the boundary into a distance r : 

b 


r 



_8_ 

21 


M 


cos <56 


(28) 


where r Q is the undistorted equatorial distance and r Q < < r fe . In Table 3 
we show the value of A r at noon ( 4 > = 0) obtained by moving the boundary from 
50 to 40 Rj due to a sudden change in the solar wind. This would increase the 
contribution of the external field near Jupiter by a factor of about 2. It is clear 
that particles adhering to a field line can easily move past Ganymede or Callisto. 
However, a sudden increase of this magnitude would move particles in by less 
than a satellite diameter at Europa, and much less than this at Io or Amalthea. 
Thus these three satellites would still act as a major barrier to diffusion. 


Electric field diffusion . Falthammar (1965), Cornwall (1968), Birmingham 
(1969), and others have stressed the possible importance of time-varying electric 
potential fields in diffusing trapped particles in towards the Earth. Falthammar 
(1968) showed that the diffusion coefficient is given by 
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(29) 


D_ 


< P n< L ’ ^,= l/r d 

2 B 2 l—i A 


where P is the power spectrum of the n th Fourier coefficient (in azimuth angle) 

n 

of the fluctuating electric field, evaluated at the drift frequency. To the extent 
that this power spectrum is independent of L, the diffusion coefficient is pro- 
portional to L 6 , instead of the L 10 dependence characteristic of magnetic 
fluctuations. 


Since it is extremely difficult, even at Earth, to estimate the magnitude of 
the electric field diffusion coefficient, we will make no attempt to do so at 
Jupiter. However, in order to evaluate the effect of a possible L 6 dependence, 
we have listed in Table 3 and plotted in Figure 5 the characteristic diffusion 
times near each of the inner satellites for a diffusion coefficient arbitrarily 
normalized to the same value as our estimated magnetic diffusion coefficient 
at L = 10, i.e., 

= 3.4 x 10" 11 L 6 R 2 /day (30) 

Such a coefficient would, of course, produce slower diffusion at L >10 but more 
rapid diffusion at L < 10. For this assumed value, diffusion times near Amalthea, 
Io, and Europa are still much longer than the calculated absorption times. 
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EFFECT OF THE DIPOLE TILT 


The calculations above of the sweeping effect of Jupiter’s satellites assumed 
that the magnetic dipole was aligned with the rotational axis. In this section we 
consider the effect of the 10-degree tilt of the dipole and calculate the resulting 
increased absorption times for near -equatorial particles. Particles which 
mirror at latitudes less than 10° (equatorial pitch angles greater than 69°), if 
located at the proper longitude, could not impact a satellite. In the co-rotating 
jovimagnetic frame of reference shown in Figure 3, a particle located at a 
jovimagnetic longitude of 0° or 180° (with the jovimagnetic prime meridian 
passing through the jovigraphic and jovimagnetic poles) would see the satellites 
passing at 10° above or below the magnetic equator. However, protons would 
drift east and electrons west (Jupiter's magnetic field is opposite in direction 
to the Earth's) with the drift periods given in Table 2. As they drift, the 
satellites appear to pass by at lower latitudes, thus interacting with some of the 
particles with equatorial pitch angles greater than 69°. At the nodes (4> m = 90° 
or 270°), even 90-degree pitch-angle particles can impact the satellites. The 
average absorption lifetime is thus a function of pitch angle. 

We first calculate the fraction of a drift period over which a particle with 
a given pitch angle can interact with a satellite. From Figure 3, the apparent 
magnetic latitude \ s at which a satellite appears to pass by can be expressed 
in terms of the magnetic longitude of the particle <p m : 

\ * 10° cos 0 

s m 
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The relationship in a dipole field between mirror latitude k m and equatorial 
pitch angle a is given by 

sin 2 a = cos 6 k /(I + 3 sin 2 \ ) 1/2 (32) 

( Hess , 1968, Eq. 2.16). For near-equatorial particles this reduces to 

a e « 90 ± (9/2) 1/2 X m = 90 ± 2.1 (33) 

Particles can interact with a satellite only if their mirror latitude k is 

m 

greater than the satellite latitude k s , or, equivalently, considering pitch angles 
between 0° and 90° only, 


(90 -<* )/2.1 > 10° (cos 4 > 


(34) 


Near the 90° node, this condition is satisfied for longitudes such that 


1 0 - 90| < sin -1 


’90 - a f 

21 


(35) 


and thus the fraction of a drift period, g , over which a particle can interact 
with a satellite is given by 


g = sin' 


'90 - a { 
, 21 


90, a e > 69° 


(36) 


g = 1, 


a e ^69° 


This fraction depends only upon the equatorial pitch angle, and is given in Table 4 
for particles of various pitch angles. (We do not consider here the case of pitch 
angles extremely close to 90°, where the finite diameter of the satellite places 
a lower limit on g.) 
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TABLE 4 


Effect of the 10-degree tilt on equatorial pitch- angle distributions. The param- 
eter g is the fraction of a drift period intheco-rotating jovimagnetic coordinate 
system over which a particle can interact with a satellite; r' bs is the modified 
absorption lifetime near the satellite Europa, as given by Equation 40; and F is 
the fraction of particles which survive absorption in diffusing past Europa, as- 
suming a characteristic diffusion time, based on magnetic fluctuations, of 24 
days. No pitch-angle scattering is assumed. 


Pitch 

Angle 

g 

Protons 

r : bs f 

Electrons 
r' F 

abs 

< 69° 

1.00 

1.9 d 

0.0016 

0.47 d 

1.2E-6 

O 

O 

c- 

0.80 

3.5 

0.011 

2.5 

0.0043 

75° 

0.51 

6.4 

0.042 

5.8 

0.034 

00 

o 

o 

0.32 

9.7 

0.086 

8.3 

0.066 

00 

cn 

0 

0.15 

17.0 

0.185 

11.5 

0.111 

o 

00 

0.09 

25.8 

0.284 

14.2 

0.147 

89° 

0.03 

67.9 

0.558 

25.1 

0.278 
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Now from Tables 2 and 3 it is clear that the absorption lifetimes for protons 
and electrons are much less than the drift periods in the co-rotating system (i.e., 
T abs < < T d)> exce P t f° r protons near Amalthea. An approximate value for the 
modified absorption lifetime r' bs can be given in each of three limiting cases: 


Case A : a e £ 69°, g ~ 1. In this case the analysis in the earlier section is 
applicable, and 


abs 


= T 


abs 


Cas e A 


(37) 


Case B: 69° < a e < 90°, but gr d >> r abs . In this case the maximum time 
that a particle would be unable to impact a satellite, assuming it were injected 
at just the right longitude, would be (1 - g) r d /2. On the average, however, a 
particle would survive about half this long (assuming, of course, that r gbs < < r d ), 
and therefore 


r abs = (1 - g) V 4 Case B 


(38) 


Case C: a = 90°, g r < < r .In this case the particles have a good 

e d abs 

chance of drifting past the nodes without impacting the satellite, and the modified 
absorption lifetime is given by 


r :bs =T abs/g Case, C 


(39) 


A semi-empirical expression for r' bs that gives nearly the correct result 
for each of these cases, yet makes a smooth transition between them, is 


T lbs = (! - g) V 4 + T abs/g AU Cases (40) 

This expression is probably reasonably accurate, within a factor of two or three, 
for all values of g. 
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In Table 4 the values of r' bs and F (as defined in an earlier section, based 
on magnetic diffusion) are given for particles of various pitch angles near the 
satellite Europa. The parameter F is plotted in Figure 7. The extremely sharp 
increase in this parameter at a e =69° is due to the fact that absorption times 
suddenly increase from those given by Equation 37 (Case A) to times which are 
of the order of a fraction of a drift period T d (Case B). 

The net effect would be a large change in the equatorial pitch-angle distribu- 
tions on either side of the orbit of Europa. On the inner side the particles would 
have distributions sharply peaked near 90°. A spacecraft which passed by the 
planet at magnetic latitudes greater than 10° would never see most of these par- 
ticles. An important parameter in the analysis of the energetic particle measure- 
ments from Pioneer 10, which will fly by Jupiter in December, 1973, will be the 
jovimagnetic latitude \ of the spacecraft. If \ jm > 10° as it moves in past 
L = 9.5, and if the processes we assume here are dominant, we would predict a 
large sudden decrease in the measured omnidirectional flux. If \ jm is near 0°, we 
would expect reduced fluxes and/or strongly-peaked equatorial pitch-angle 
distributions at L < 9.5. 

At Io and Amalthea, however, even though absorption lifetimes for particles 
with pitch angles greater than 69° would be substantially increased, these life- 
times are still very much shorter than the characteristic diffusion times given 
in Table 3. Thus, if the basic assumptions of this study are correct, and if the 
J calculated diffusion times are reasonably accurate, even near- equatorial particles 
cannot diffuse past Io or Amalthea. 

However, we know from the analysis of the decimeter radiation that electrons 
with sharply-peaked pitch-angle distributions are present near L = 2 ( Thorne , 1965), 
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AT EUROPA r djff = 24 DAYS 


PROTONS 



ELECTRONS 


EQUATORIAL PITCH ANGLE (a e ) 


Figure 7. The fraction _F of particles that can diffuse in past Europa without being absorbed, 
plotted as a function of pitch angle. The characteristic diffusion lifetime is assumed to be 
24 days. The increase at 69° < a e < 1 1 1 0 is due to the fact that during a substantial fract 
of these particles drift period, the satellites appear to pass by at latitudes above or below 
their mirror points. 





just inside the orbit of Amalthea. If there exists some kind of mechanism which 
diffuses electrons much more rapidly than we have calculated, those with pitch 
angles near 90° would have a greater probability of surviving than those with 

pitch angles less than 69°, since the drift period for 50 Mev/gauss electrons 
near Amalthea is about 17 days, compared with an absorption lifetime (Case 2) 
of 2.4 days. This effect could easily produce the kind of pitch-angle distributions 
required by Thome. 


45 



DISCUSSION 


A number of additional simplifying assumptions were made in deriving the 
results of the earlier sections. Here we examine the consequences of removing 
some of these simplifications. 

The existence of quadrupole and higher-order terms or a possible offset 
of the equivalent dipole would not affect our analysis in any significant way. 

The effect of a dipole offset would be to make the circular satellite orbits appear 
to be offset by an equivalent amount in the co-rotating jovimagnetic frame of 
reference shown in Figure 3. However, each satellite would come back to the 
same apparent position after one co-rotation period. For random absorption 
(Case 1), the same particles would be subject to absorption again (assuming no 
diffusion or longitudinal drift), and for snowplow absorption (Case 2), exactly 
the same path would be plowed out during each apparent orbit. A similar argu- 
ment applies if we consider the presence of quadrupole or higher-order terms. 

We can also argue that the slight eccentricity in the orbit of each satellite 
is unimportant. From the eccentricities listed in Table 1 we calculate that 
the maximum deviations from a circular orbit are 540, 0, 200, and 1600 km, 
respectively, for satellites Amalthea through Ganymede. For all but Amalthea 
these are much less than a satellite diameter. 

The alternate compression and extension of a co-rotating field line by the 
solar wind during one jovian day is another possible source of radial motion 
of particles. Equation 28, giving the intersection of field lines with the magnetic 
equator, can be used to calculate the maximum diurnal deviations at noon ( 4 > = 0) 
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and midnight ( 4 > = 180°) from their average positions. These distances are 0.3, 

22, 220, and 2280 km, for satellites Amalthea through Ganymede, respectively, 
assuming r b = 50 Rj. Again, these distances are much less than the respective 
satellite diameters, and can be ignored. 

We have assumed that all trapped particles co-rotate with the planet. If at 
some radial distance this assumption does not hold, the apparent period of revolu- 
tion of each satellite would be the sidereal period given in Table 1, rather than 
the co-rotation period. This would mean that the absorption of 50 Mev/gauss protons 
would be governed by Case 2 instead of Case 1 at Europa and Ganymede, and the 
absorption times for both electrons and protons at Io, Europa, and Ganymede 
would increase to 1.8, 3.6, and 7.2 days, respectively, instead of the values 
given in Table 3. This relatively small increase would not have a significant 
effect on our findings. 

We have assumed in this study that the source of Jupiter's radiation belts 
are solar wind particles which have diffused in from the outer magnetosphere. 

We have ignored possible additional source mechanisms, such as CRAND 
(cosmic ray albedo neutron decay). The flux of galactic cosmic rays incident 
on Jupiter's atmosphere is reduced by a factor of approximately 4000 from the 
corresponding flux at Earth, since the cosmic-ray cutoff ridigity at Jupiter is 
360 times as large. Only cosmic-ray particles with energies greater than about 
5 Tev (5 x 10 12 ev) can reach the surface of Jupiter at low latitudes. Thus the 
CRAND source strength will be several orders of magnitude less than at Earth. 

The loss rates will also be much lower, due to the lower density of thermal 
plasma and essentially complete absence of neutral atmosphere above a thousand 
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kilometers. The resulting steady-state flux of trapped particles from CRAND 
will probably be no greater than, and probably less than, the flux of CRAND- 
produced particles at Earth. As with inward-diffusing particles, Amalthea at 
L = 2.5 and Io at L = 6 will act as barriers to the outward diffusion of these 
particles, unless our calculated diffusion times are substantially in error. 
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SUMMARY AND CONCLUSIONS 


In this study we have examined the competing processes of (1) absorption of 
trapped particles by impact with one of Jupiter’s satellites and (2) diffusion of these 
particles by violation of the third adiabatic invariant. We have assumed that a 
particle would be removed from the' radiation belt if, in its normal motion, it 
would happen to impact a satellite. We have assumed that the satellites are un- 
magnetized and that they would not interfere with the motion of a nearby particle 
in any way. With the additional assumption that Jupiter’s dipole axis is aligned 
with its rotational axis, the average absorption lifetime of a proton or electron 
with magnetic moment ft = 50 Mev/gauss was found to be 0.5 - 5 days near 
Amalthea, Io, Europa, or Ganymede, and somewhat longer at Callisto (see Table 3). 

A diffusion coefficient was then calculated, assuming violation of the third 
invariant due to magnetic fluctuations associated with fluctuations in the solar 
wind. We assumed that on the average there was one disturbance per (Earth) 
day which doubled the contribution of the external field near Jupiter. This led to 
a diffusion coefficient D equal to 0.13 r 10 /r® day -1 = 3.4 x 10" 15 L 10 R^/day, 
taking the subsolar boundary distance r b = 50 Rj. A similar calculation for the 
Earth's magnetosphere, taking r fe = 10 R E , gives D = 1.3 x io -9 L 10 R^/day, 
about equal to the value quoted by Williams (1972) as that value which could 
explain the behavior of steady-state and non- steady-state outer-belt protons 
at Earth. 

A characteristic diffusion time in the vicinity of each satellite, equal to the 
mean time for a particle to diffuse by one satellite diameter , was computed. 


49 



This ranged from 137 years at Amalthea to 0.002 days at Callisto (see Table 3 
and Figure 5). This time was short compared to an average lifetime at Callisto, 
roughly equal to the absorption time at Ganymede, but substantially greater than 
the absorption time at Europa, Io, and Amalthea. The fraction of particles which 
survive absorption in diffusing past each satellite was computed (see Table 3). 
Amalthea and Io were found to act as complete barriers to diffusion, and Europa 
would be expected to reduce particle fluxes by three orders of magnitude or more. 
A small reduction would be observed at Ganymede, and Callisto would have 
essentially no effect. Figure 6 illustrates the predicted proton fluxes at Jupiter, 
showing the expected decreases at L = 6 and L = 9.5. Similar calculations were 
made for a diffusion coefficient proportional to L 6 which might be produced by 
fluctuating electric fields, arbitrarily normalized to the value of the magnetic 
diffusion coefficient at L = 10. The conclusions were essentially unchanged. 

We then examined the effect of the 10-degree tilt of Jupiter’s magnetic 
dipole axis. The absorption times were found to depend strongly upon equatorial 
pitch angle. For particles with pitch angles greater than 69°, mirroring at lati- 
tudes less than 10°, the average absorption times are increased from those 
calculated earlier to times which are a fraction of a drift period (see Table 4), 
since at some jovimagnetic longitudes near-equatorial particles cannot interact 
with any satellite. Figure 7 shows the fraction of particles that would survive 
absorption in diffusing past Europa as a function of pitch angle, assuming a 
characteristic diffusion time at Europa of 24 days. The equatorial pitch-angle 
distributions of particles just inside this satellite would thus be expected to be 
highly peaked near a e =90°. 
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However, even for near-equatorial particles, the absorption lifetimes near 
Io and Amalthea are substantially shorter than diffusion times based on magnetic 
fluctuations produced by solar wind variations. Thus, given the assumptions of 
this study, solar wind particles diffusing inwards from the outer magnetosphere 
are not expected to be able to penetrate the inner satellites of Jupiter, particularly 
Amalthea and Io. Yet high fluxes of trapped electrons with energies of the order 
of 10 Mev are definitely responsible for the synchrotron emission near L = 2. Are 
these solar-wind electrons, and if so, how did they get there? There are several 
possible explanations: 

1. A satellite might in some way perturb the motion of a nearby electron, 
so that it would not impact the satellite. This could be due to an inherent mag- 
netic field possessed by the satellite, or to the interaction of a highly-conducting 
satellite with the field lines which co- rotate past it. The strong effect which Io 

is observed to exert on the decametric radiation bursts indicate that it interacts 
strongly with the field lines at L = 6, giving some support to the latter hypothesis. 
Protons with the same magnetic moment as electrons will have a relativistic 
momentum everywhere 43 times greater (the square root of the mass ratio), and 
therefore the motion of protons is not likely to be affected nearly so much by 
local perturbations in the magnetic or electric fields. It is possible, therefore, 
that electrons could slip past the inner satellites but protons would be absorbed. 

2. Diffusion of electrons in the inner magnetosphere might be much more 
rapid than the calculations based on magnetic fluctuations have indicated. If 
Amalthea and Io do not perturb the motion of nearby electrons, diffusion rates 
in this region would have to be increased by about four orders of magnitude for 
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electrons to diffuse past these satellites without being absorbed. Additional 
processes, unknown at Earth, might be responsible for this rapid diffusion at 
Jupiter. Brice (1972) and Brice and McDonough (1972) have suggested that 
neutral atmospheric and ionospheric winds and turbulence at the top of Jupiter's 
atmosphere might provide the required driving mechanism. 

If indeed the diffusion of electrons proceeds much more rapidly than indicated 
in this study near Io and Amalthea, the effect of the 10 -degree tilt of the dipole 
will be for these satellites to preferentially remove particles with mirror lati- 
tudes greater than 10 degrees. Those that survive will have pitch-angle distribu- 
tions strongly peaked near 90°; we believe that this is the most probable mechanism 
producing the sharply-peaked distributions required in the studies of Thorne 
(1965) and Clarke (1970), particularly since Amalthea at L = 2.5 is located just 
outside the peak of the synchrotron emission region. Without more precise 
knowledge of the processes causing enhanced diffusion, however, we cannot say 
whether protons would be affected by the same rapid diffusion. 

3. Additional unknown acceleration processes might give rise to energetic 
electrons in the inner magnetosphere. We know that energetic electrons are some- 
times injected at low L-values at Earth following a major magnetic storm (e.g., 
Williams et al., 1968). Similar processes might operate at Jupiter. It is not as 
likely, however, that energetic trapped protons could be produced by such means. 

We probably will not be able to determine which alternative best explains 
the existence of energetic electrons at L = 2, or whether substantial fluxes of 
protons are able to diffuse past the satellites, until actual particle measurements 
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are made at Jupiter by Pioneer 10 and later spacecraft. It will be very important 
for the onboard detectors to distinguish clearly between trapped electrons and 
protons, as there is good reason to believe that they might behave very differently. 
It will be particularly important to get accurate flux measurements near L= 5.95, 
9.47, and 15.1, the positions of Io, Europa, and Ganymede. Accurately-measured 
pitch-angle distributions, particularly if the flyby is at low magnetic latitudes, 
will be extremely useful in evaluating the various hypotheses. 

Finally, if a Jupiter flyby is planned into the very inner regions of its 
magnetosphere, and if for the safety of the spacecraft the lowest possible radia- 
tion levels are desired, the mission should be planned so that the spacecraft 
remains at relatively high magnetic latitudes, so as to avoid intense fluxes of 
particles with equatorial pitch angles near 90°. 
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